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Linear water wave theory suggests that wave patterns caused by a steadily moving 
disturbance are contained within a wedge whose half-angle depends on the depth- 
based Fronde number Fh- For the problem of flow past an axisymmetric pressure 
distribution in a finite-depth channel, we report on the apparent angle of the wake, 
which is the angle of maximum peaks. For moderately deep channels, the depen¬ 
dence of the apparent wake angle on the Fronde number is very different to the 
wedge angle, and varies smoothly as Fh passes through the critical value Fh = 1. 
For shallow water, the two angles tend to follow each other more closely, which 
leads to very large apparent wake angles for certain regimes. 


The linear dispersion relation for small amplitude water waves gives rise to intricate 
and distinctive surface wave patterns that have been the subject of enduring interest to 
applied mathematicians and physicists since the 19th century^. When these patterns form 
at the wake of a steadily moving object, they are often referred to as ship wave patterns 
or ship wakes, regardless of whether or not the object is an actual ship. Relatively simple 
geometric argument^ show that ship wave patterns propagate within a travelling wedge; 
for an infinitely deep fluid, this wedge has a half-angle of ^wedge = arcsin(l/^ ^ 19.47° 
(referred to here as the Kelvin angle) regardless of the speed of the disturbance. 

For a fluid of finite depth, the derivation of the wedge angle is more complicatecP, and 
turns out to depend on the speed of the disturbance U, which itself is measured by the depth- 
based Fronde number Fh = U/\/gH. Here H is the depth of the fluid and g is acceleration 
due to gravity. Of particular interest here, the wedge angle has the following properties: 
^wedge ^ arcsin(l/3)~^ as Fh ^ ^wedge ^ Tij‘2 as Fh ^ 1 •> ^wedge ^ Tij2 as 
Fh -t 1+; ^wedge ^ as Fh -t oc. We note that for Fh < 1, the flow is referred to as 
being subcritical. Subcritical flows are characterised by two types of waves, transverse and 
divergent. Transverse waves run roughly perpendicular to the direction the disturbance is 
travelling. Divergent waves propagate at an acute angle to this direction. For supercritical 
flows, Fh > 1, the wave pattern is entirely made up of divergent waves. 

There has been a resurgence of interest in ship wave patternJ^, with a particular focus on 
the wake angle that is actually observed in practic^. Here we pay attention to what we call 
the apparent wake angle, 6>app, which is the angle between the line passing through the high¬ 
est peaks of the wave and the centreline. For axially symmetric pressure distributions acting 
on the free surface of an infinitely deep fluid, numerical simulations show that the apparent 
wake angle is slightly less than the Kelvin angle arcsin(l/3) for moderately fast distur¬ 
bances; however, for sufficiently fast-moving disturbances, the apparent wake angle appears 
to decrease like the inverse of the length-based Fronde number, Fl = U/\/^ L, where L is a 
length-scale associated with the pressure distributiorP^ (see also Ref^^^^. Further recent 
studies have been undertaken on the effect of non-axisymmetric pressure distrib utionj^dJlMJ^ 
a continuous distribution of sourceJI^, consta nt vor ticit^fl^, nonlinearit}ff^^I^ and surface 
tensiorP^. A comparison of the results in Ref j^ l ^F^ ig given by He et alW^. All of these 
studies are for an infinite-depth fluid. 

By considering the problem of flow past an axially symmetric pressure distributions in 
a fluid of finite depth, we shed light on the seemingly contradictory results for the wedge 
angle and recently reported apparent wake angle. As noted above, the wedge angle ^wedge 
increases from arcsin(l/3) to 7r/2 as Fh increases to Fh = 1, and then monotonically 
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decreases for Fh > 1- However, for a moderately deep fluid, we find the dependence of 
the apparent angle 6>app on Fh is very different, and in fact behaves much like the infinite 
depth case. This is remarkable because the transition between subcritical {Fh < 1) and 
supercritical {Fh > 1) hows is characterised by distinct changes in qualitative behaviour 
in the wave pattern and the singularity structure of the dispersion relationship, while the 
apparent wake angle 6>app changes smoothly as Fh increases through Fh = 1. On the other 
hand, for shallow water, the apparent angle ^app follows the wedge angle ^wedge closely, 
giving rise to very large apparent wake angles for moderately small length-based Froude 
numbers, Fl. 

In the following, we consider an axisymmetric pressure distribution of strength Pq and 
intrinsic length L acting on the surface of a fluid flowing in the positive x direction with 
speed U. We will present solutions for a fluid of infinite depth and a finite depth H. 
To enable a direct comparison with recent worli^, our dimensional pressure distribution is 
p(x, y) = Pq exp(—For the infinite depth problem we choose to scale speeds 
by U and lengths by L. This leads to two nondimensional parameters, €l = Po/{pgL), a 
measure of pressure strength (where p is the fluid density), and the length-based Froude 
number, Fl. Using linear theory, an exact solution for the free-surface C(x, y) in the infinite 
depth case isPIl 


C{x,y) = - CLp{x,y) + ^ J j 

— 7r/2 0 


7r/2 oo 

Cl f f k‘^p{k^'ip) cos{k[\x\cos'ip F ysin^p]) 


k — k() 


dk dp) 


eLH{x) 
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— 7r/2 


where p{k^ p) = ^^exp(—(47r^))/7r is the Fourier transform of our pressure distribution 
and H{x) is the Heaviside function. The path of /^-integration is diverted below the pole 
k = ko, where ko = 1/(P| cos^ p). 

For the finite depth problem we instead use H as the length scale, which leads to three 
nondimensional parameters: the pressure strength €h = Po/{pgH), the depth-based Froude 
number Fh, and the ratio of pressure length scale to fluid depth S = L/H. The nondimen¬ 
sional parameters for the finite and infinite depth problems are related by Fl = Fnj^ 
and tL = The exact solution for the finite-depth problem i^^ 


= - ^HP{x,y) + 
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where po = 0 for Fh < 1 and po = arccos(l/Fi^) for Fh > 1. As with the infinite depth 
solution Q the path of /^-integration is diverted below the pole, which this time is denoted 
hy k = kpp), where ki is the real positive root of 

kPn — sec^ p tanh k = d, pQ < p (3) 

In both the infinite-depth 0 and finite-depth 0 solutions, the double integral rapidly 
tends to zero far away from the pressure distribution. Therefore, to analyse the wake in the 
far field, only the single integral term is required. 

To measure the apparent wake angle, 6>app, we need to locate the line that passes through 
the highest peaks of the wave, and then calculate the angle between that line and the 
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centreline ^ = 0. The method for locating the highest peaks depends on the flow regime. 
For subcritical flows (and the infinite depth case), we isolate each transverse wave, mark 
the highest point, and fit a line throu gh th e highest points of all the wavelengths (we call 
this method 1, used previously in Refs^®, for example). For supercritical flows, there are 
no transverse waves, so we instead use the highest points on each of the divergent waves to 
fit the line (method 2). 

In Fig. we present data for apparent wake angle against both Fronde numbers Fl and 
Fh for a variety of dimensionless pressure length scales 5. These data are represented by 
the (red) solid dots. As just mentioned, the data for the infinite depth case and for the 
subcritical flows Fh < 1 are calculated with method 1, while data for supercritical flows 
Fh > 1 are calculated with method 2. We see in Fig. [^b)-(d) that for moderately deep 
channels, the qualitative behaviour is very similar to the infinite depth case in Fig. dla) 
considered previousl}!^ that is, the apparent wake angle 6>app is approximately arcsin(l/3) 
for Fl less than some value (roughly Fl = ^ or slightly less), then decreases like F^^ for 
large values of Also plotted in each panel of Fig. is ^wedge (a dashed line). It is 
interesting to see that for deep to moderately deep channels, the apparent wake angle 6>app 
does not follow 6>wedge- 

It is insightful to observe the wave patterns themselves in Fig. especially as the flow 
transitions from subcritical Fh < 1 to supercritical Fh > 1- For sufficiently deep channels 
{S sufficiently small), the transition has almost no effect on 6>app (also refer back to Fig.f^b)- 
(d), where the slight kink in the data for at Fh = 1 is due to the change in method 
of measuring the angle). In Fig. [^a)-(b), we present wave patterns for a typical small 
value (taken to be S = 1/4), computed for Fh slightly less than and slightly greater than 
the critical value Fh = 1. The corresponding plan view and centreline plots are shown in 
panels (e)-(f). In the plan view, the highest peaks are represented by open circles, while the 
angle ^wedge is represented by the dashed lines. These surface profiles help to show how, for 
sufficiently deep channels, the apparent wake angle 6>app varies smoothly through Fh = 1 
and is not at all close to 6>wedge- 

Returning to Fig. for shallow channels, as in Fig. [^e)-(f), there are clear differences 
between supercritical and subcritical flows. In the former regime {Fh > 1, ^ 1), the trend 

in the data is dominated by the F^^ decay. While for the latter regime {Fh < 1, (5 ^ 1), 
the apparent wake angle 6>app can be very large (much larger than arcsin(l/3)), and here 
^app does follow ^wedge closely. The transition between subcritical and supercritical flows 
has a dramatic qualitative effect on the free surface (Fig. |^c),(d)(g),(h)) with the absence 
of the transverse waves for Fh > 1, exemplified in the centreline plots of Fig. |^g),(h). 

To calculate the asymptotes for the wake angle O^pp in Fig. (the thick blue line and 
the thin green line) we perform a far-field stationary phase approximation on the second 
integral in ([^ using polar coordinates (r, ^) (defined via x = rcos^, y = rsin^) to obtain 

C(r, 0) « a(r, 0) sin(r^(7/i(6>), 0) - 7r/4), 
where ^(t/, 0) = ki{%jj) cos('0 + 0) and 

n{r ft) = exp (-(5^fci(V>i(6>))V(47r^)) 

(-F'ff-sec(?/)i(6»))2sech(A:i(?/)i(6»)))2)y5^;^;(^07(^)7^ 

Recall ki is the real positive root of (§. The function ipi{0) satisfies 

IT 

k[{'tpi) cos(V’i + 6>) - fci(V’i) sin(V’i + 0) = 0, tpo < tpi < - and 0 < 0 < ^wedge, 

while is the second partial derivative of g with respect to t/. 

We must now find the maximum of a{r^0). The two approaches depend on the method 
used to measure the wake angle ^app- For method 1, r is held constant (to r = 1, say) and 
we then compute the turning point of a(l,^). For method 2 we must choose r such that 
^hi{rg{%jji{0)^0) — 7r/4) = 1, which leads to the new amplitude function 

const exp (-(5^fci(V>i(6>))^/(47r^)) / 

(F^ - sec{'tpi{9))Hec\i{ki{'tpi{0))Y y 9,p,p{^i{9),9)' 


tV’) 


( 5 ) 
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FIG. 1: Measured apparent wake angles in degrees (solid circles) plotted against the 
Fronde number for the linear solutions Q and for a fluid of (a) infinite depth and 
(b)-(f) finite depth, presented on a log-log scale. The dashed line is (a) the Kelvin angle 
arcsin(l/3) ~ 19.47° and (b)-(f) the wedge angle ^wedge in the finite-depth case. The thick 
(blue) line is the theoretical asymptote given by method 1 (valid for subcritical flows) 
and the thin (green) line is asymptote given by method 2 (required for supercritical flows). 


A simple numerical test shows that, for each value of S, the maximum of both 0 and 
^ behaves like const as Fh oc, so that, in both cases. 




Fh 




VSFl 


as Fh oc. 


(6) 


The numerical values of P for both methods are found to increase from ^ = 0 at (5 = 0 to 
P ^ 1~ as(5^oc, thus the asymptotic behaviour of the wake angle is dominated by finite 
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(c) Fh = 0.95, 5 = 9 {d) Fh = 1.5, 5 = 9 



(e) Fh = 0.95, 5 = 1/4 (f) Fh = 1.05, 5 = 1/4 


(g) Fh = 0.95, 5 = 9 


(h) Fh = 1.5, 5 = 9 


FIG. 2: (a)-(d) Free surface profiles for various values of Fh and 5. (e)-(li) A plan view of 
the profiles in (a)-(d) with a centreline plot. The open circles are the highest peaks and 
the solid line is the line of best fit. The dashed lines the wedge angle ^wedge- 


depth effects and approaches the finite depth wedge angle 6>wedge- As shown in Fig. we 
note that method 1 gives [3/^/6 l/(40^/^v^) as (5 ^ 0, which agrees with the infinite 

depth result computed by Darmon et al.^^. 

Thes e resu lts are worth comparing with the phenomenological model of Moisy and 
Rabaud^^^, which for the present finite-depth problem can be extended in a straightfor¬ 
ward fashion to predict that 

tan6>app = c^(l - - c^Cp) with /c = 27r/(5, (7) 

where Cg and Cp are the group and phase velocities, respectively, defined by Cg = doo/dk 
and Cp = uj/k^ and the dispersion relation is cd = {k tanh k)^^/F h- As Cg = 0{F^^) and 
Cp = 0{F^^) as Fh oc, we see that Moisy and Rabaud’s method gives ^app ^ Cg with 
k = 27 t/6 which, using the notation in (|^, gives 

(5tanh(27r/^)-1-27rsech^(27r/(5) 

2y^27r6 tanh(27r/(5) 

This result is also shown in Fig. (dashed line). In the deep water limit, ^ gives 
(2v^)“^ as (5 ^ 0, which agrees with ReP^. For shallow water ^ ^ 1, it collapses to our 
method 2. 

We now plot in Fig. |^the apparent wake angle ^app against Fh for fixed values of Fl- 
Note that for each fixed value of Fl, an increase in Fh must be associated with an increase in 
6 so that Fl = Fh / V^. If we imagine a given pressure distribution with dimensional length 
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FIG. 3: Plots of (a) [3 and (b) [3/\/^ in ^ against S. The thick (blue) line conies from 
method 1, the thin (green) line from method 2, the dashed (purple) line by 


scale L moving at dimensional speed I/, then each set of data points in Fig. [^corresponds 
to decreasing the channel depth from infinite {Fh = 0) to zero {Fh = oo). With this in 
mind, FigureQa) shows that for a faster moving stream (see data for Fl = 1.5 and 2.5), the 
transition from subcritical to supercritical flow via a decreasing fluid depth has little effect 
on the apparent wake angle ^app- Ultimately, regardless of the speed of the stream (that is, 
for all fixed values of Fl)^ the apparent wake angle asymptotes to the wedge angle ^wedge 
for sufficiently shallow depths (see Figure j^b) for large Fh)^ It is interesting to compare 
these results with the phenomenological model §; the agreement is not so good for deep 
channels, but improves as the depth decreases. 




FIG. 4: Measured apparent wake angles 6>app in degrees plotted against the depth-based 
Fronde number Fh for the linear solution <§ for constant values of Fl = 0.7 (blue 
squares), 1.5 (red dots) and 2.5 (green triangles). Note that as Fh increases, 5 increases 
such that Fh /= Fl . The solid lines represent Q , while the dashed line indicates the 

wedge angle 6>wedge- 


In summary, we have considered the prototype problem of flow past an axisymmetric 
pressure distribution acting on the surface of a finite-depth channel. For the linearised 
solution, we have explored the dependence of the apparent wake angle 6>app (or angle of 
maximum peaks) on the dimensionless variables. This work extends previous studies for flow 
past an axisymmetric pressure distribution for the infinite depth case. It has been insightful 
to fix the depth of the channel and increase the speed of the flow (fixed 6^ increasing Fh^ see 
Figure and also to fix the speed of the flow and decrease the depth (fixed increasing 
Fh^ see Figure 1^. It is remarkable that for much of the parameter space, the apparent 
wake angle ^app varies smoothly as Fh is increased through critical value Fh = 1. On other 
hand, for shallow water, we see the apparent wake angle 6>app tends to follow the wedge 
angle 6>wedge closely, which for near critical flows produces very large apparent angles, much 
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higher than the Kelvin angle arcsin(l/3). 

We close by noting that the we have deliberately used a very simple pressure distribution, 
which more easily facilitates a comparison with many recent studies, but may not accurately 
reflect the wave patterns produced by a real ship. An alternative simplified approach, for 
example, is to employ a pair of point pressures on the surfacd^ (one positive and one 
negative), or a continuous distribution of sourced. Indeed, the wave pattern generated by 
a finite-depth flow past a pair of point pressures was considered very recently by Zhu et 
aijlSl study found that for subcritical flows in very shallow channels, the wake angle 

was close to the wedge angle 6>wedge: which is very similar to our observations. On the 
other hand, their argument based on interference effects suggests that the wake angle scales 
like in the large Fronde number limit, which is very different to (^. This difference 
in the large Fronde number scaling is, at least in part, likely due to th e difference in the 
flow configurations considered (see the relevant discussions in Ref j^^ l ^^ * about the effect of 
non-axisymmetric pressure disturbances on the large Fronde number scaling). Either way, 
a further advantage of our approach here is that we provide direct comparisons between 
analytical predictions and measured values. 
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